Abstract. The main scope of this paper is to investi- 
gate the possible existence of a metallicity dependence 
of the overshooting from main sequence stars turbulent 
cores. We focus on objects with masses in the range 
~ 2.5Mq— ~ 25Mq. Basically, evolutionary time scale 
ratios are compared with star numbers ratios on the main 
sequence. Star populations are synthesized using grids 
of evolutionary tracks computed with various overshoot- 
ing amounts. Observational material is provided by the 
large and homogeneous photometric database of OGLE 
\ 2 project for the Magellanic clouds. Attention is paid to 
• the study of uncertainties: distance modulus, intergalactic 
| and interstellar reddening, IMF slope and average bina- 
^ ■ rity rate. Rotation and chemical composition gradient are 
£h ' also considered. The result for the overshooting distance is 
i I™? = 0.40i°;J 2 6 H p (Z - 0.004) and l™° = 0.10t° ;^H p 
(Z = 0.008) suggesting a possible dependence of the ex- 
tent of the mixed central regions with metallicity within 
the considered mass range. Unfortunately it is not yet pos- 
■ sible to fully disentangle effects of mass and chemical com- 
OO ' position. 
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1. Introduction 

Large convective phenomena occur in the cores of main se- 
quence stars with masses above about 1.2 Mq (for galactic 
chemical composition). In standard models, convection is 
crudely modeled with the well-known Mixing Length The- 
ory of Bohm-Vitense (1958) ) (hereafter MLT) and the core 
extension is determined according to the Schwarzschild 
criterion. The Schwarzschild limit is the value of the ra- 
dius where the buoyancy force vanishes. However inertia 
of the convective elements leads to an extra mixing above 
the Schwarzschild limit, called "overshooting" and usually 
expressed in fraction of the pressure scale height. Several 
theoretical works (for a review see Zahn 1991) give argu- 



ments in favor of such an additional mixing. Many labo- 
ratory experime nts show evidences for overshooting (see 
Massaguer 1990| ). Although overs hooting can occur below 
an external convective zone (see Alongi et al. 1991 ), this 
paper is exclusively concerned with core overshooting. 
One of the first empirical determination of convective 



core overshooting was obtained by Maeder & Mermilliod 
|(1981| ) who used a set of 34 galactic open clusters and fit- 
ted the main sequence width with an additional mixing 
of about 20-40 % in mass fraction. Mermilliod & Maeder 



(1986 ) derived an overshooting amount of about 0.3 H p 



again for solar- like chemical composition and for a 9-15 
M range. |5tothers fe Chin (199l"l ) derived an overshoot- 
ing amount < 0.2 H p for Pop. I stars using the metal- 



enrich ed opacity tables published in Rogers fc Iglesias| 
1992| ). 



During the last decade, many evolutionary model grids 
have been computed with an overshooting am ou nt equal 
or close to 0. 2 H p : e.g. Charbonncl ct al. (1996 ) or Bcrtelli| 
ct al. (1994). This second team uses a formalism (see 



Brcssan et al. 1981) slightly different from the Geneva 



team one. Generally the same overshooting amount is used 
whatever the metallicity and mass are. 

Kozhurina-Platais et al. (1997) obtained l OVCI ~ 0.2 ± 
0.05 Hp for the galactic cluster NGC 3680 (solar metallic- 
ity) with the isochronc technique, this method consists in 



fitting the cluster CMD features ( particularly the turn-of f 
position) with model isochrones. |lwamoto fc Saio (1999 ) 



compared evolutionary models with observations of three 
binary systems: V2291 Oph, a Aur and 7/ And ("binary 
system" technique). The authors adjusted either the he- 
lium content or the overshooting parameter to get a better 
fit to observations. The best results were obtained with a 
moderate overshooting amoun t (^ 0.15 Hp). For su per- 
solar metallicity (Z Q = 0.024) [Lebreton ct al. (200l| ) de- 
rived l OVCI < 0.2 Hp from the modeling of the Hyades 
cluster turn-off. 



Maeder & Mermilliod (1981) have suggested an over- 
shooting increasing with mass within the studied range of 
2-6 Mq which is also found by Schroder et al. (1997 ) with 
a study of binary systems. According to their results, the 
overshooting should increase from 0.24 Hp for 2. 5 Mq 
to ^ 0.32 H p for 6.5 Mq. With a similar study |Ribas 
|et al. (2000| ) also found a mass dependence. 

The question of a metallicity dependence must also be 
addressed. Ribas et al. (200C )'s results suggest a slight 
metallicity dependence for a stellar mass around 2.40 Mq 
(see their Table 1). The more metal poor star SZ Cen 
(in mass fraction: Z = 0.007) is satisfactorily modelled 
with an overshooting distance 0.1 H p < l ov < 0.2 Hp 
and objects with Zq ranging between 0.015 and 0.020 



seem to have an overshooting around 0.2 Hp. Keller et al 
|(200l| ) have recently explored the dependence of over- 
shooting with metallicity by means of the isochrone tech- 
nique using isochrone grids from the Padova group. Their 
study involves HST observations of four clusters: NGC 
330 (S MC), 1818, 2004 and 2001 (LMC). |Keller et al 
|(2001| ) find the best fit (with respect to age and over- 
shooting) for an overshooting amount which is equiva- 
lent to l OVCT = 0.31 ± 0.11 
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In this paper, we carry out an independent study of 
a possible metallicity dependence of overshooting with a 



technique which differs from the "binary system" (Riba; 
[et al. 200C ; Andersen 1991) and "isochrone" techniques 
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Our method is based on star-count ratios, with compar- 
isons between observational material and synthetic pop- 
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ulation results in color-magnitude (CMD) diagrams. We 
are then led to discuss several points: particularly distance 
modulus, reddening and binarity rate. If the dependence 
of overshooting with metallicity (or mass) was thereby to 
be firmly assessed, it would then be a challenge to under- 
stand its physical origin. 

We are concerned with a metallicity range relevant to 
the Magellanic Clouds and take advantage of the homo- 
geneous OGLE 2 data, which provide color magnitude di- 
agrams for ~2x 10 6 stars in the Small Magellanic Cloud 
(hereafter SMC) and ~7x 10 6 in the Large Magellanic 
Cloud (hereafter LMC). On the theoretical side, we es- 
timate the number of stars from evolutionary model se- 
quences computed with different amounts of overshooting. 
From these data sets and using evolutionary models with 
intermediate and low metallicity, we estimate the over- 
shooting value during the main sequence in the SMC and 
LMC for a stellar mass in the range 2.5 Mq - 25 Mq. 

In Sect. |we describe the observational data involved 
in this work. Sect. |3| is devoted to the method used: data 
selection and star counting. Sect. || gives the main fea- 
tures of our population synthesis procedure. Sect. is de- 
voted to astrophysical inputs, and Sect. || to results and 
effects of uncertainties. Sect. discusses the results. It 
must be emphasized that we determine in fact the extent 
of the inner mixed core region which can be due either to 
true overshooting or to another process as rotation, some 
observational evidences exist about correlation betweenn 
metallicity and vs'mi (see Venn et al. 1999), the prob- 
lem of rotation is shortly discussed in the Sect. [?]. Finally, 
Sect. |§| gives some comments and concluding remarks. An 
appendix has been added to provide details about the pop- 
ulation synthesis algorithm and error simulations. 

2. Observational data 

The observational data set considered here has been ob- 
tained by the Optical Gravitational Lensing Experiment 
(OGLE hereafter) consortium during its second operat- 
ing phase (for more details and references the reader can 
consult URL: jittp : //www, astrouw. edu.pl/~ogle/ ). 



2.1. SMC and LMC data 



We have downloaded the SMC data described in Udalski 



et al. (1998). The data used in this paper are from the 
post-Apr. 8, 2000 revision. The SMC is divided into 11 
fields (labeled SCI to SC11) covering 55'xl4'; each field 
contains between ~ 100,000 and ~ 350,000 objects. For 
each object several quantities are available: equatorial co- 
ordinates, BVI photometry and associated standard errors 
<tb, (t y and 07. This database has the great advantage of 
being extensive and very homogeneous. 

The LMC data are described in |Udalski et al. (2000| ). 
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The BVI map of the LMC is composed of 26 fields (SCI to 
SC26) in the central bar of the LMC. The dataset includes 



Fig. 2. (a) Data from SMC with a v < 0.015 mag and 
box definitions (Ni + N 2 = 4653 and N 2 /N 1 = 1.08), 
(b) the same for the LMC case (Ni + N 2 = 4113 and 
N 2 /N 1 = 1.01). 



photometry and astrometry for about 7 million stars over 
a 5.7 square degree field. 



3. The star-count method 

3.1. Data selection 

As shown in Fig. ^b, the standard error on V-magnitude, 
cry, increases with the magnitude. This is also true for B 
or I-magnitudes. Hence the errors on (B — V) or (V — I) 
colors rapidly increase and reach values as large as 0.2 
mag around a V-mag^ 20: this is of the same order as the 
Main Sequence width. 

As we are interested in the MS structure and as we 
must minimize error effects while keeping a quite good 
statistics, we have chosen to take into account only data 
with ay and ob (or a\) lower or equal to 0.015 mag, lead- 
ing to a maximum error on color of 0.02 mag. The value 
of 0.015 mag appears to be an optimal choice maintaining 
a good statistics with photometric errors remaining small 
compared with the MS width. Fig. [l] sketches the pro- 
posed selection process and displays differences between 
the entire Color-Magnitude Diagram (Fig. |l|a) and the fi- 
nal diagram (Fig. ^c): obviously, the remaining data are 
those corresponding to lower magnitudes. 

This selection process leaves ~ 4700 objects on the 
SMC MS (over a total of more than 2.2 millions objects) 
in the BV system (~1100 objects in the VI system) and ~ 
4000 objects on the LMC MS (over a total of more than 
7.2 millions objects) in the BV system (~1600 objects in 
the VI system). As we can see the BV system presents a 
more favorable statistics, therefore in the following we will 
work only with this set of bands. 

Tables 4 from |Udalski et al. (1998^ |2000| ) indicate that 
completeness for V < 18 should be better than about 99 % 
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Fig. 1. (a) CM-Diagram for the OGLE 2 SC 1 field, 10 % of the data have been plotted for sake of clarity, (b) Standard 
deviation of the measurements: ay versus magnitude V, the dot-dashed line indicates the limit- value considered (0.015). 
(c) The resulting CM-diagram after selection (all the fields within SMC have been plotted). 



for the SMC; and should be around 96 % - 99 % depending 
on the field crowding for the LMC. 

3.2. Star count ratios: an observational constraint 

As the absolute number of stars arriving on the ZAMS 
per unit of time for a given mass is unknown, we rather 
compute star count ratios. To count stars, we first define 
an area in the CM-diagram. As we are interested in the 
MS structure, we choose a region which contains the main 
sequence "bulge" revealed after the data selection process 
(see Fig. ||a) with the most convenient geometrical shape: 
a "parallelogram" (for automatic count purpose). A couple 
of opposite sides (AB and DE in Fig. ||a) are chosen to be 
more or less "parallel" to the main sequence axis. 

In the CM-Diagram, main sequence stars evolve from 
the blue to the red side. The MS width is mainly an evolu- 
tionary effect connected to a characteristic time scale tms 
(time spent by a star on the Main Sequence). The distri- 
bution of the objects within the Main Sequence should be 
related to this time scale. Therefore we divide our paral- 
lelogram into two regions called "box 1" and "box 2" (see 
Fig. ||a) where the respective numbers of objects N\ and 
N 2 are similar (N 2 /Ni ~ 1). This ratio is taken as an ob- 
servational constraint and it will enable us to discriminate 



between theoretical grids of evolutionary tracks computed 
with various overshooting amounts. 

We now turn to the method used to build a synthetic 
stellar sample comparable to the OGLE 2 ones (after se- 
lection) from evolution simulation outputs. 

4. Population synthesis 

4-.1. Evolutionary models 

Our evolutionary model s are built with the ID Henyey 
type code CESAMQ (see |Morel 1997) ) in which we brought 
several improvements. Applying modern techniques like 
the projection of the solutions on B-spline basis and au- 
tomatic mesh refinements, CESAM allows robust, stable 
and highly accurate calculations. We use as physical in- 
puts: 



the OPAL 96 opacities from Iglesias fc Rogers (1996[) at 
high temperatures (T > 10000 K) and the Alexander 



& Ferguson (1994) opacities for cooler domains. For 



metallicity higher than the solar one (that occurs dur- 
ing the He core burning phase) we use elemental opac- 
ities (Los Alamos) calculated by Magee et al. (1995). 



1 CESAM : Code d'Evolution Stellaire Adaptatif et Modu- 
laire 
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the EFF equation of state from Eggleton et al. (1973) 
elemental abundances are from Grevesse & Noels 
(19[H (the "GN93" mixt ure), the cosmological helium 
is from [zotov et al. (1997j ): Yp = 0.243, and the helium 



content is scaled on the solar one following a standard 
helium-metallicity relation: Y = Y P + Z(AY/AZ). 
The calibration of a solar model in luminosity yields 
AY/AZ = 2 (|Lcbrcton et al. 1999|) from the calibra- 



tion of the solar model radius. This value is com patible 
wit h the recent v alue AY/AZ = 2.17 ± 0.40 of |Pcim 
ber |t et al. (2000|) . We therefore adopt AY/AZ 



2. 



ftp : //tangerine . astro .mat .uc . pt/pub/BaSeL/. This 
library provides color-calibrated theoretical flux distribu- 
tions for a large range of fundamental stellar parameters, 
T cS (2000 to 50,000 K), log 5 (-1.0 to 5.5 dex), and [Fe/H] 
(-5.0 to +1.0 dex). The B&SeL flux distributions are cali- 
brated on the stellar UBVRIJHKL colors, using: 

— empirical photometric calibrations for solar metallicity 

— semi- empirical relations constructed from the color dif- 
ferences predicted by stellar model atmospheres for 
non-solar metallicities. 



— for the chemical composition we adopt [Fe/H] der ived 



fro m Ccpheid measurements by Luck et al . (1998 i 



Details about the calibration procedure are given in Leje 



for the SMC they find a range from -0,84 to 



q g^, unc et al. (1997 ) and Lejcune et al. (1998). Compared to 



with a mean value: [Fe/H] = —0.68 which leads to 
X Q = 0.745, Y„ = 0.251 and Z = 0.004. 
— for the LMC they find a range from —0.55 to —0.19, 
combining all the values we obtain a mean value of 
[Fe/H] = -0.34 leading to: X = 0.733, Y = 0.259 
and Z = 0.008. 
The nuclear reaction rates are from |Caughlan fc 
Fo oler (19881), except: 12 C(a,7) 16 Q, 17 0(p, 7 ) 18 F 
taken from jCaughlan et al. (1985| ) and 17 0(p,a) 14 N 
taken from Landre et al. (1990). The adopted rate for 
12 C(a,7) 16 is quite similar to the NACRE compila- 
tion ( Angulo et al. 1999]) one: a factor of about two 



higher than Caughlan & Fowler (1988) and about 80% 
of Caughlan et al.'s (1985) one. 

To take into account the metallicity effect on the mass 
loss rate ( |de Jager et al. 19*88 ) we adopt the scaling fac- 
tor (Z p /0.02) u tl derived from the Kudritzki fc Hummei| 
(19|^) models. 

The convective flux is computed according to the clas- 
sical MLT. We use a mixing length value Imlt = 
1.6 Hp. Th is 
et M- (1992j) 



value has been derived by Schaller 
from the average location of the red 
giant branch of more than 75 clusters. Very similar 
val ue (1.64) ha s been found more recently by Lebre- 
ton et al. (1999 ). An extra-mixing zone is added above 
the Schwarzschild convective core: this "extra-mixing" 
zone is set to extend over the distance Z ovcr = a ver Hp , 
clover being a free parameter, the value of which is dis- 
cussed here. 

— the external boundary conditions are determined in a 
layer within a simple grey model atmosphere built with 
an Eddington's T(r) law. 

4-2. Conversion of the theoretical quantities into 
observational ones 

In order to compare theoretical results to obser- 
vational data, conversions are needed. Transforma- 
tions of the theoretical quantities, (Mboi, T c g) into 
absolute magnitudes and colors are derived from 
the most recent version of the Basel Stellar Li- 
brary (B&SeL, version 2.2), available electronically at 



the previous versions of the B&SeL library, all the model 
spectra of stars with T e ff > 10, 000 K are now calibrated 
on empirical colors from the T e g versus (B — V) relation 
of Flower (1996). In addition, the calibration procedure 
for the cool giant model spectra has been extended in the 
present models to the parameter ranges 2500 K < T c ff < 
6000 K and -1.0 < \ogg < 3.5 f\ 

4-3. Population synthesis 

In contrast with "classical" works on population synthesis 
where the CMD as a whole is simulated, we construct a 
small part of the CMD: the area containing the brighter 
MS stars. In this way the task is simplified. Artificial stel- 
lar samples have been generated from our evolutionary 
tracks with a specially designed population synthesis code 
CReSyPS|. 

In our framework the main hypothesis is that the Star 
Formation Rate (SFR) is constant during the time scales 
involved here: i.e. a few hundred megayears. So for a given 
mass the number of observed stars (i.e. those correspond- 
ing to a given evolutionary track) must be proportional to 
the time scale of the main sequence. We assume that the 
SFR is constant in time and mass (equal for all masses in 
the range explored in this work), if we note r the SFR: 
At w 1/r represents the mean time elapsed between two 
consecutive star births. For the observational star samples, 
At is unknown but the objects numbers are available. We 
choose At to get similar total star numbers in boxes 1 and 
2 (i.e. iVi + N2) both in the synthetic CMD and observa- 
tional diagram. We point out that the ratios N2/N1 are 
not sensitive to the At value chosen. 

The evolutionary track grids scan a mass range be- 
tween 2.5 M and 25 M from the ZAMS to logT off ~ 3.8 
covering the entire box ranges in color and magnitude (de- 
fined in Sect. |3.2| ) . The mass step is increasing from 0.5 
Mq around 3 Mq stars to 5 Mq above 15 Mq. Sev- 
eral overshooting amounts have been used from 0.0 to 



2 In the previous versions of the B&SeL models, we adopted 
T ff = 5000 K and log o = 2.5 as the upp er limits for the 



calibration of giants (see Lejeune et al. 199; 



Code Rennais de Synthese de Populations Stellaires 
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Fig. 3. Synthetic CM-Diagrams for SMC and LMC chemical compositions, panels (a) and (b) are for the SMC with 
two overshooting amounts: a over = 0.0 and a over = 0.8, panels (c) and (d) are for the LMC with: a over = 0.0 and 
clover — 0.2 respectively. For panel (a) only 30% of the synthetic objects have been displayed for clarity purpose 
and corresponding evolutionary tracks have been plotted. For other panels: (b), (c) and (d) the number of displayed 
objects has not be reduced. In all cases the total number of stars in the boxes -used for calculations- is close to the 
observational one: (a) Nx + N 2 = 4403 {N 2 /N 1 = 0.65), (b) JVi + N 2 = 5011 {N 2 /N 1 = 1.59), (c) ATi + N 2 = 3582 
(N 2 /Ni = 0.94), (d) m +N 2 = 4717 {N 2 /N! = 1.13). We recall that empirically we got for the SMC Nx +N 2 = 4653 
{N 2 /N! = 1.08) and for the LMC Ni + N 2 = 4113 (N 2 /N 1 = 1.01). The ratios N 2 /Nx obtained theoretically are rather 
indcpcndant from the value of N\ + N 2 in the synthetic CM-Diagram, for instance in the case of the panel (b) we 
got N 2 /N 1 = 1.60 with N x + N 2 = 6399. The cloud of dots in panel (b) (l over = 0.8 H p , SMC) calls a comment: it 
appears to be bimodal, i.e. showing over-populated regions around V <~ 17.5 mag and V ~ 16.3 mag. Indeed for high 
overshooting values the main sequence of masses as small as ~ 3Mq can reach V ~ 17.5 generating with their large 
evolutionary time scale an over-populated region. Moreover the binarity shifts a part of this population to a ~ 0.8 
mag brighter region, with a binarity rate (3 = 0.0, this "bimodal effect" disappears 
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0.8. CReSyPS treats the photometric errors by simulat- 
ing OGLE 2 ones (see App. |A|) which is very important 
for our purpose. Our algorithm requires the knowledge of 
some input parameters: distance modulus, reddening and 
absorption, binarity rate, Initial Mass Function (hereafter 
IMF) slope and photometric errors. 

We summarize here the main steps of the algorithm: 



5.1. Distance modulus 

Large Magellanic Cloud. The LMC distance modulus has 
a key role in extragalactical distance determinations, but 
its value is still debated. The de terminations range be- 
tween "short" distance scal es (i.e. Btanck et al. 1998 ) and 
"long" distance scales (i.e. Laney fc Stobie 1994). Using 



the HIPPARCOS calibrated red clump stars, Btanek et al 



ST |EP 1: a mass distribution is generated between 2.5 (1998|) found ^q.lmc — 18.065 ± 0.031 ± 0.09 mag and 



Mq and 25 Mq following the Salpeter's law: dN/dm 
m -a Sa i P (- see Sect. 5.4). 

— STEP 2: for each mass, an evolutionary track is inter- 
polated within the grid calculated by the evolutionary 
code. On each track, models are selected every time 
step At, which is adjusted in order to yield a total 
number of stars equivalent to the observed one. 

— STEP 3: consistently with the value of the binary 
rate < (3 > (see Sect, p 7 ^ ), objects are randomly se- 
lected to belong to a binary system and the magni- 
tudes of these systems are calculated. Triple systems 
(and higher multiplicity systems) are neglected. 

— STEP 4: distance modulus is added (and in the case of 
SMC a random "depth" inside the cloud) and synthetic 
photometric errors are attributed to magnitudes (see 
App. 

— STEP 5: we use a "quality filter": objects with too 
large photometric errors are rejected from the syn- 
thetic sample. 

— STEP 6: color is calculated, reddening and extinction 
coefficient are applied. Concerning reddening, a gaus- 
sian distribution is applied around the mean value in 
order to simulate object-to-object variations (see dis- 
cussion in Sect. |5.2| ). 

With the interpolation between evolutionary tracks, it 
is very important (particularly at low mass, i.e. 3.0Mq < 
M < 4.0Mq) to reproduce the time scale tms with a good 
accuracy. A test at 3.25 Mq has shown that the "interpo- 
lated time scale", T^ ipo1 , is very close to the calculated 
one (with the evolutionary code) t^I with a difference 
not larger than about 1%. Also important are the magni- 
tude interpolations on the Main Sequence: our tests also 
show a very good agreement between interpolated mag- 
nitudes and calculated ones, differences are unsignificant 
(about 10~ 3 — 10~ 2 mag, whereas the photometric errors 
are much larger). 

Our code intensively uses a random number generator. 
We have chosen an algorithm insuring a very large period 
about 2 x 10 18 (program "ran2" from Press ct al. 1992Q , 
which is much larger than the number of synthetized ob- 
jects. 

As a result, examples of synthetic samples generated 
by CReSyPS are displayed in Fig. || where the influence of 
overshooting is shown for both clouds. 

5. Astrophysical inputs 



Laney fc Stobie (1994 ) from a study of Cepheids Period- 
Luminosity relation obtained /iq,lmc — 18.53 mag with an 



internal error of 0.04 mag. Grocncwcgen fc Salaris (200l[ ) 
found ^o,HV2274 = 18.46 ± 0.06 mag from a study of the 
LMC-eclipsing binary system HV 2274. They indicate a 
LMC center distance at ^o.lmc = 18. 42 ± 0.07 mag. Re- 
cently, from the DENIS survey data, Cioni et al. (2000) ) 
derived a distance modulus for the LMC of /*o ; lmc = 
18.55 ± 0.04 (formal) ±0.08 (systematic) using a method 
based on the apparent magnitude of the tip of the red giant 
branch. The HST Key Project Team adopted /l*o,lmc = 
18.50 ± 0.15 mag ( [Mould et al. 2000| ). In order to bracket 
the most recent estimations, we have chosen the HST Key 
Project value: 

Mlmc = 18.50 ± 0.15 mag 



Van der Marel fc Cioni (200l| ) give an order of magnitude 
of the depth of the LMC. They indicate small corrections 
to magnitude for well studied individual objects within the 
LMC, ranging between A^ ,lmc = -0.013 (SN 1987A) to 
A^o,lmc = +0.015 (HV 2274). We neglect these correc- 
tions which have the same order of magnitude than the 
photometric errors. 

Small Magellanic Cloud. Laney & Stobie (1994) suggest 
a distance modulus (based on Cepheids) of /xo.lmc = 
18.94 mag with an internal error of 0.04 mag; this mod- 
ulus decreases by about 0.04 mag if ca librators are half- 
weighted. More recently Kovacs (2000) (with a method 
based on dou ble mode Cepheids ), find ^ ,lmc = 19.05 ± 
0.13 mag and Cioni et al. (2000| ) have ^o.lmc = 18.99 ± 
0.03 (formal) ±0.08 (systematic) mag. We retain the fol- 
lowing estimation: 

/iSMC = 18.99 ±0.10 mag 

The SMC distance modulus only represents an average 
distance, prowl et al. (2001 ) have evaluated the depth of 
the SMC along the line-of-sight by a study of populous 
clusters. They derived a depth between ~ 6 kpc and ~ 
12 kpc; these values lead to magnitude differences of 0.2 
and 0.4 mag respectiv ely. Previous studies, see for instance 
pardiner et al. (1991 ) , show similar results with a line-of- 
sight SMC depth ranging between ^4 — 7 kpc and ~ 15 
kpc strongly depending on the location in the SMC. 

We have chosen to model the SMC depth with a gaus- 
sian distribution of distances around ^smc with a stan- 
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dard deviation: 



Fig. 4. Relative number of stars N st ars/N to t versus red- 



a 



depth 

SMC — 



0.05 mag 



kpc (about ~ 0.3 



which represents a total depth of 
mag). 

5.2. Reddening and absorption 

We have to distinguish: foreground reddening 



E(B — IV )mw (clue to material m Milky Way) and 
internal reddening ]^(_h5 — V); with an origin into the 
Cloud itself. These quantities are expected to change 



dening from Dutra et al. (2001) (dashed curve), OGLE 2 
experiment Cepheid catalogue: see |Udalski et al. (1999a| ) 
and Udalski et al. (1999b), and data read in Fig. 24 of 
pestreicher ct al. (1995[ ) (dot-dashed curve), (a) Data for 
SMC, (b) data for LMC. 

E(B - V) MW +i = 0.15 ± 0.02 (LMC). In Fig. | we have 
displayed the histogram of E(B — V)mw+i values from Du 



tra et al. (2001) and OGLE group. OGLE data have a 



along tltic hnc-ol-sight. Here we model the total reddening 
as E(B - V) M w+i = E(B - V) M w + ETB - V) ; taking 
into account its non-uniformity. From the literature, we 
derive estimations for the mean value and the dispersion 
of E(B — V)MW+ii object-to-object variations can then 
be simulated. 

We now discuss reddening determinations for SMC and 
LMC. From a study of spectral prop erties of galactic nu- 
clei behind the Magellanic Clouds, Dutra et al. (2001 ) 
have evaluated the foreground and background redden- 



better statistics with respectively 1333 (SMC) and 2049 
(LMC) obje cts, against 14 (SMC) and 22 (LMC) for |Dutra 
ct al. (200ip . Dutra et al's data are systematically less red, 



this could be inherent to their method: they observed ob- 
jects behind Clouds and observations are easier through 
the more transpar ent regions of the clouds. 



In addition, Ocstrcichcr et al. (1995| ) have de- 



termined the reddening for 1503 LMC foreground 
stars with a UBV photometry based method: 
E(B — V)mw = 0.06 ± 0.02 mag, a quite low value 
because it is related to foreground stars. It shows a 



ings fot both Clouds. For the LMC, they found an average 
spectroscopic reddening ot ETB — V)mw+i = 0.12 ± 0.10 
mag. The uncertainties essentially come from the de- 
termination of the stellar populations belonging to 



spread (0.02) similar to the OGLE 2 one. Pestreicher 



background galaxies: in the case of LMC, when Du- 



tra et al. (2001| ) consider only red population galax- 
ies, they find E(B — V)mw+i = 0.15 ± 0.11 mag, which 
gives an idea of the global uncertainty on E(B — V), 
which should be around ~ 0.02 — 0.03 mag (about 
- 13 - 20%). For the SMC putra et al. (200l| ) find 
E(B - V) M w+i = 0.05 ± 0.05 mag. 



et al. (1995 ) distribution is in very good agreement (see 
Fig. |4|)) with Dutra et al's one which tends to confirm 
that Dutra et al's result could be underestimated (Dutra 
et al's results are supposed to take account for foreground 
and internal reddeing). Therefore in the case of LMC, we 
prefer to retain the OGLE average value for purpose of 
consistency: 

< E(B - V)^ c +i >= 0.15 mag 

Fig. |i|b shows that the distribution shape is the same for 
Dutra et al. (200l|) and |Oestreicher et al. (1995|) , OGLE 
one being quite narrow which appears slightly underesti- 
mated, thus we take a value similar to Dutra et al. one: 



0.5 - 
0.4 - 
0.3 - 
0.2 - 
0.1 - 
{ 



SMC Data 1 



(a) 



T LMC 

E(B— V)mw4 
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(b) 
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We take into account an additional uncertainty on < 
E(B - V)LMC + . > of about: 

<5e(b-v) = 0-02 mag 

The SMC case is more questionable (we only have two sets 
of data), we favor the OGLE values because they are likely 
more suitable for performing simulations which synthesize 
OGLE data. Moreover OGLE data have larger statistics. 
We adopt: 

< E(B - Vjg^+i >= 0.09 mag 
with a crudely estimated uncertainty of about: 
<*I(B- V) = 0-015 mag 

In this case also, the OGLE standard deviation (see 
Fig. p|a) seems to be low, therefore we adopt the Dutra 



CcpheM — ste* — diocovorod — m — both — Cloudo. — OGLE 

values are: E(B - V) M w+i = 0.09 ± 0.01 (SMC) and 



et al. (2001) one: 



7 SMC 
E(B— V)mw4 



0.05 mag 
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The absorption coefficient is taken from Schlegel et al 
(19980 

A v = 3.24 x E(B - V) M w+i 
and is calculated for each object. 



5.3. Binary rate 

Evaluating the average binary rate < (3 > in objects as ex- 
tended as the Magellanic Clouds is not easy. Locally (i.e. 
within a particular area of the galaxy) this multiplicity 
rate depends -at least- on two factors: (1) the star den- 
sity and the kinematics of the objects which influence the 
encounter probability; (2) the initial binary rate (relative 
number of binaries on the ZAMS) . Within the Magellanic 
Clouds, the binary rate likely varies over a wide range and 
we only consider its spatial average value < (3 > . 

Chez (1995) finds in solar neighbourhood that for main 
sequence stars and young stars the binary rate < /3 > 
ranges between 0.10 and 0.50 (it peaks at < /3 >= 0.50). 
Therefore we tested the effects of binarity for these two 
extreme values. 

In our population synthesis code, binaries are taken 
into account with a uniform probability for the mass ratio 
q = A/2 /Mi (in the considered mass range). 

5.4. IMF slope 

The IMF has been extensively discussed by many authors. 
Toward both Galactic poles and within a distance of 5.2 
pc from the Sun, Kroupa et al. (1993) found a mass func- 



tion: dN/dm « m QSalp with asai p ~ 2.7 for stars more 



massive than 1 Mq. In the LMC, |Holtzman et al. (1997| ) 
inferred -from HST observations- a value consistent with 
the Salpeter (1955|) one: asai p ~ 2.35. At very low metal- 
licity, Grillmair et al. (1998| ) observed the Draco Dwarf 
spheroidal Galaxy ([Fe/H] « -2) with the HST. They 
concluded that the Salpeter IMF slope remains valid in 
the Magellanic Clouds and we have chosen: 



Q!Salp 



2.30 ±0.30 



However we must keep in mind that some circularity in 
work exists when using an IMF. As described by Garcia 



fc Mermilliod (200l| ) the IMF can be derived from the 
observed Present Day Mass Function (PDMF) using evo- 
lutionary tracks and their corresponding time scales which 
depend on the adopted value for the overshooting ! 

5.5. Star Formation Rate 

For a given mass, the Star Formation Rate (SFR) rep- 
resents the number of stars "created" per unit of time. 
Vallenari et al. (1996| ) have studied three stellar fields of 
the LMC and have found a time scale of about ^2—^4 
Gyr for the "bulk of star formation" . We therefore make 
the reasonable assumption that the SFR remained quite 



constant during the short galactic period relevant for this 
work, i.e. for the last ~ 300 Myr. The SFR involved here 
is an average value over each cloud. 

6. Resulting overshooting amounts 

6.1. Large Magellanic Cloud 

As a first step we choose the mean values for each astro- 
physical input (discussed in Sect. [5]), this yields for the 
LMC the following overshooting: 

Uer = 0.09 Hp 

which is a rather mild amount. We examine in Fig. ^[ 
how the Z 01 , er -value is affected by the uncertainties on the 
astrophysical inputs: 

— Changing the IMF slope asaip in the range 2.0 — 2.6 
we obtain: 

0.02 < S < 0.09 Hp 

which tends to minimize the overshooting. 

— Next, a test with the average binary rate < (3 > in the 
range 0.10 — 0.50 leads to: 



0.00 < 1<P> < 0.14 Hp 

— A distance modulus value in the range 18.35 < /xlmc ^ 
18.65 enables us to derive: l OV er — 0.0 Hp (in fact for 
A*lmc = 18.35 all the values for simulated N2/N1 are 
larger than the observed one) and I over = 0.21 Hp for 
Mlmc = 18.65. 

— An average reddening between 0.13 and 0.17 leads re- 
spectively to lover — 0.0 Hp (in this case also all the 
values for simulated N^/Nx are larger than the ob- 
served one) and l OV er = 0.27 Hp. 

We stress that uncertainties on distance modulus and 
reddening infer the largest uncertainties on the final over- 
shooting values. We retain for the LMC average chemical 
composition: 



l LMC = 10 

'over — u.J-U 



+0.17 



^-o.io H P 

which indicates that a mild overshooting amount around 
~ 0.1 — 0.2 Hp is needed to model LMC stars as found 
in the majority of determinations involving solar chemical 
composition objects (see Sect. |l|). 

6.2. Small Magellanic Cloud 

For the SMC, using the mean value of each astrophysical 
inputs we obtain (see also Fig. ||): 

Z OV er = 0.41 Hp 

— If the IMF slope varies between extreme values 
(—2.0 < asaip ^ ~2.6), the overshooting varies within 
the following boundaries: 



0.36 < &E. <0.50H 



P 
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Fig. 5. Overshooting determinations for SMC (panels a, b, c, and d) and LMC (panels e, f, g and h). The influence 
of distance modulus, reddening and IMF slope are considered for each cloud: continuous lines correspond to central 
values of these parameters discussed in Sect. |E| and dashed lines to the associated error bars. Inferred values of i ver 
(and uncertainties) are given on each figure. For (c) and (g) panels, dotted line is for < (3 >= 0.0 and dot-dashed 
for< j3 >= 0.5 (see text). 



similarly, an average binary rate ranging between 0.10 
and 0.50 leads to: 

0.34 < l<£> < 0.45 Hp 

— The uncertainty on SMC distance modulus (18.89 < 
Msmc ^ 19.09) leads to: 

0-37 < Cer < 0.46 Hp 

— Similarly if one considers the uncertainty on the aver- 
age reddening (< E(B — V) > ranging between 0.075 
and 0.105), the overshooting amount shows a high sen- 
sitivity to reddening: 

0-36 < lSS -v) < 0.52 Hp 

Again, uncertainties on distance modulus and reddening 
are the largest. We retain for the SMC: 

-,+0.12 



;SMC 
^over 



0.40t£& Hp 



Whatever the simulation is, statistical errors are of the 
order of 0.01 Hp which can be safely neglected. In the 
SMC case, the required overshooting appears to be much 
larger than for LMC stars and for solar composition stars. 



7. Discussion 



7.1. An upper limit with Roxburgh's criterion 



Roxburgh's criterion ( Roxburgh 1989] ) is a very general 
constraint on the size of the convective core. It is written 
as an integral formulation over the stellar core radius: 



-R co 



(,-^rad ^nuc] . dr — 



r=0 



1 dT 



47rr 2 c!r(l) 



where L ra d and L nuc are respectively the radiative en- 
ergy flux and the total energy flux (in J.s -1 ) generated 
by nuclear processes, r is the radius, R cor e is the core size 
including the "overshooting" region. $ represents the vis- 
cous dissipation (in J.s _1 .m -3 ). In the whole stellar con- 
vective core the turbulence is supposed to be statistically 
stationary and the temperature gradient has to be almost 
adiabatic. In Eq. ^ the integrand is positive when r is lower 
than the Schwarzschild boundary where L rac i — L nuc and 
it becomes negative beyond. 

The viscous dissipation $ is unknown but the integral 
constraint is satisfied for larger R cor e value when $ = 0. 
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Table 1. Time weighted average overshooting distances 
for a 6 Mq main sequence model, derived with the 
"Roxburgh's criterion" neglecting dissipative phenomenon 
($ = 0). 



Metallicity Z 
Average EOA 



0.004 (SMC) 0.008 (LMC) 



0.6 Hp 



0.6 H D 



Hence , iieglecLing the dissipation by selling <f> — piu- 

vides the maximum possible extent of the convective core 
which can be considered as the upper limit for overshoot- 
ing. Evolutionary tracks have been calculated, using Rox- 
burgh's criterion, for a representative mass of 6 M Q and 
SMC and LMC metallicities. The equivalent overshooting 
amount (EOA), given in Table 1, is the time weighted av- 
erage overshooting distance along the evolutionary tracks, 
expressed in pressure scale height. 

In both cases (LMC and SMC), Roxburgh's criterion 
predicts a maximum value (i.e. neglecting viscous dissi- 
pation) around 0.6 Hp (see Tab. |l|) independent from 
Z . Our determinations -i.e. Z™ r c = 0.40t°;Jg H p and 



LMC 



0.10 



+0.17 
-0.10 



Hp- therefore are compatible with the 
theoretical upper limit given by the Roxburgh's criterion. 



7.2. Influence of rotation 

In addition to convection, rotation is an other important 
phenomenon inducing mixing through shear effects and 



other instabilities. For instance Venn (199£ ) finds surface 
abundance variations in SMC A supergiants that could be 
explained by some kind of mixing related to rotation. 

Taking account of the rotational effect brings new im- 
portant unknown features: (1) the f2-value distribution 
and (2) the vsmi distribution for the considered stellar 
population. Both features remain unconstrained by obser- 
vational studies. 

In addition, stellar rotation involves many effects and 
physical processes that are non-trivial to include in mod- 



ern evolutionary codes. Talon et al. (1997) show that (see 



their Fig. 5) a rotating ID-model with an initial surface 
velocity of 300 km.s -1 leads to a main sequence track 
equivalent to an overshooting model using l oveI = 0.2 Hp. 
Despite great theoretical eff orts, a free parame ter remains 
for horizontal diffusivity in Talon et al. (1997 ) treatment 
of rotational mixing (see Zahn 1992 ). 

Rotation changes the global shape of an evolutionary 
track, through two distinct effects: (1) the material mix- 
ing inside the inner part of the star which brings more 
fuel into the nuclear burning zones like overshooting, (2) 
the effective surface gravity modification leading to color 
and magnitude changes (which depend on the angle be- 
tween the line-of-sight and the rotational axis) . In their 
Fig. 6, Maeder fc Meynet (200l| ) show the influence of ro- 
tation on evolutionary tracks for low metallicity objects 



(Z = 0.004). These tracks have been calculated taking 
into account: (1) an "average effect" on surface, (2) the 
internal mixing. These tracks are very similar to those 
calculated with different overshooting amounts values. 

An additional effect which needs to be discussed here is 
the surface effect: modifications of colors and magnitudes 
of MS stars due to rotation (in absence of any mixing 



phenomenon) have been studied by Maeder & Peytremann 
(1970D with uniformly rotating models. Their Tab. 2 gives 
expected changes of My and (B — V) as a function of 
(angular velocity expressed in break-up velocity unit) and 
vsmi (this latter ranges from to 457 km.s -1 , for a 5 
Mq star). In this table standard deviation for My and 
(B — V), arc: om v = 0.18 mag and o-(b-v) = 0.01 mag. 
Therefore the rotation effect has roughly the same order of 
magnitude than present uncertainties on magnitudes and 
colors. 



Data from Wolff & Simon (1997) 
(15 <m < 18.5) 



sin i (km.s ) 



Fig. 6. Relative number of st ars (N st ars/N to t versus V x 
sin i for Wolff fe Simon (1997 ) objects with magnitudes in 
our studied range. 



However, stars with vsmi greater than ~ 200 km.s 1 are 



quite rare, as shown by the data of Wolff & Simon (1997) 
(see Fig. |6|). Then keeping only data with vsmi < 200 
km.s -1 , leads to: gm v — 0.20 mag and cr(s_y) = 0.001 
mag. The effect on absolute magnitude remains of the 
same order, whereas the effect on color becomes largely 
negligible. We conclude that our results remain valid, even 
if the major mixing is due to rotation. In this case, the 
value of Z over would change its meaning. Major contribu- 
tion to i ver value would represent a shear effect mixing. 



7.3. Influence of chemical composition gradient 

We have sofar assumed a uniform chemical composition. 
The chemical composition may vary inside each Magel- 
lanic Cloud. The existence of an abundance gradient in 
the Clouds is still debated and spectroscopic measure- 
ments with a statistics as large as the statistics of OGLE 
2 data are not available. In their Tab. 4, Luck et al 



data are not available. In their Tab. 4 
|(1998| ) give spectroscopic determinations of [Fe/H] for 7 
SMC Cepheids and 10 LMC Cepheids. For SMC data, 
the standard deviation is cp^/H] ~ 0.07 dex leading to 
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negligible variations for the heavy elements mass frac- 
tion Zq. Therefore the SMC can be considered as chemi- 



cally homogeneous for our purpose. For LMC, Luck et al 
(1998|)]frnd a standard deviation ^rp^/H] ~ 0.10 ^ex gi ym S 
0.007 < Zq < 0.01. From evolutionary tracks of typical 
mass (6 Mq) and an overshooting of 0.1 Hp, changing Z 
from 0.007 to 0.01 has a negligible effect on magnitude 
and an effect of ~ 0.003 mag on color, which is largely 
lower than the photometric errors. We conclude that -in 
the light of the present knowledge- the chemical composi- 
tion gradient does not change our results significantly. 



7.4- Comparison with other works 

From the investigation of young clusters in the Magellanic 
Clouds, Keller et al. (2001) did not find any noticeable 
overshooting dependence with metallicity. They obtained 
for NGC 330 (Z ~ 0.003) C?r C33 ° = 0-34 ± 0.10 H p , 
which is compatible with our determination for the SMC: 



Z£mc = QAO+m Hp. For NGC 2004 (Z 



0.007) Keller 



ct al. (2001| ) got C?r C2 °° 4 = 0-3 1 ± 0- 11 H p: whilc for 



similar metallicity we derived /qvcV — O.IOIq io Hp which 
is also compatible with Keller et al.'s result. One can note 
that masses involved in our simulations (average mass of 
~ 7 — 8 Mq with a standard deviation of 4 Mq) are higher 
than the Keller ct al. (200l| ) one (terminus masses in the 
range 9-12 Mq for the four clusters). [Keller et al. (2001| ) 
do not discuss the influence of the uncertainty on distance 
modulus of the clusters and use (J-lmc — 18.45 mag and 
Msmc = 18.85 mag. 

Ribas et al. (2000) derive overshooting amounts from 
evolutionary models of galactic binary systems. For SZ 
Cen (Z Q ~ 0.007) they find 0.1 < l OVCI < 0.2 H p which 
is close to our value for the LMC, but the mass of SZ 
Cen is 2.32 Mq and some mass effect cannot be avoided, 
therefore any comparison with the present results must be 
considered with care. 



; 4 


4mU - — *—& « 


1 . 1- 4 f - 




► « 1 



_ u 

In Fig. we summarize results from several authors. De- 
spite the small number of points, a slight dependence of 
overshooting with metallicity cannot be excluded. How- 
ever, at low and high metallicities, the considered mass 
ranges are different and the errors remain substantial, 
therefore a definite conclusion is not yet possible. 



Fig. 7. Overshooting parameter a ovet versus metallicity 
Zq from various sources. Open triangles represent results 
from [Ribas et al. (2000[ ) for SZ Cen (Z ~ 0.007) error 
bars have been indicated, arrows mean that the derived 
value is a minimum. The open square shows a result from 



Kozhurina-Platais et al. (1997) for the galactic cluster 



NGC 3680, error bars are indicated. Filled triangles are de- 



terminations from Keller ct al. (2001): with continuous er- 
ror bars amounts corresponding to the SMC cluster NGC 
330, NGC 1818, NGC 2100 and with dashed error bars 
result for the LMC cluster NGC 2004. Open circle: de- 
termination in Hyades cluster from Lebreton et al. (2001 ) 
(upper limit for overshooting). Filled diamonds: SMC and 
LMC determinations performed in this work. Errors on Z 
have been evaluated assuming an internal error on [Fe/H] 
ofO.ldex. 



8. Conclusion 

In this paper we have estimated the overshooting dis- 
tance from a turbulent core for intermediate-mass main se- 
quence stars. The result for SMC is Z^MC = OTOt^ Jg H P , 
and for the LMC Z^mc = O.lOT^fo h p- The main contri- 
butions to errors are those brought by distance modulus 
and reddening uncertainties. We have shown that chemi- 
cal gradients within the clouds and rotation surface effects 
of studied stars cannot significantly influence our results. 
Binary rate, IMF slope have no important effects as well. 
For SMC, despite different methods and data, we find a 
result very similar to Keller et al.'s (2001) one for cluster 
NGC 330. The case of LMC is more questionable because 
of the rather large uncertainty on reddening. 

Fig. ^ tends to indicate a sensitivity of overshooting to 
metallicity. However a mass effect cannot be excluded, we 
can only stress that if such a dependence exists, it should 
be an increase of overshooting with decreasing metallicity. 
However, the overshooting is expected to increase with 
mass, unfortunately studied samples at solar metallicity 
have often lower masses than those at low metallicities. 
Therefore further investigations arc needed to disentangle 
these effects. In any cases, if this dependence is confirmed 
the next challenge will be the physical explanation of this 
metallicity-overshooting effect. 

Finally, the overshooting amounts derived in this work 
have a statistical meaning: they are average values over 
time (in real star "overshooting" likely changes during the 
main sequence) and over mass in the considered range. 
Moreover these amounts represent an extra-mixing above 
the classical core generated either by inertial penetration 
of convective bubbles or shear phenomena related to rota- 
tion. The real extent of the core likely results from a com- 



bination of both processes, indeed rotating models Maedei 



Meynet (2001 )'s rotating models still need overshoot- 



ing. 
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Appendix A: Photometric error simulations 

As we selected the data using a criterion involving the pho- 
tometric standard deviation of magnitude measurements, 
we have to generate an artificial standard deviation for the 
theoretical magnitude computed from evolutionary mod- 
els. Moreover the general properties of the synthetic stan- 
dard deviation distribution must be similar to the OGLE 
2 one. 

We describe here the scheme used to generate the 
pseudo-synthetic photometric standard error distribu- 
tions. The prefix "pseudo" means that we have extracted 
information about the standard error distribution from 
the OGLE 2 data themselves (see Fig. |A.l| (a)). For that 
purpose, we divide the relevant range of magnitudes into 
bins; in each bin, we construct the histogram of standard 
deviation values (Fig. A.l(b)). This histogram then is fit- 
ted with a function of the form: 
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Fig. A.l. (a) Standard deviation cry versus V-magnitude 
for objects belonging to the SCI field of the SMC. (b) His- 
togram of the cry values for magnitude V between 20.5 and 
21.0, the fit (dashed curve) is performed with a function 
of a type given in the text. Differences between the fit and 
the histogram are clearly insignificant for our purpose, (c) 
Synthetic try distribution generated with our algorithm. 

cr TO value is too large) or the magnitude to is changed into 
m nmsv , following a gaussian distribution having a stan- 
dard deviation a m - 



Let us comment about differences between Fig. A. 1(a) 
and Fig. |A.l| (c). Fig. A. 1(a) contains the "evolutionary in- 
formation" -i.e. more objects at high magnitudes- whereas 
Fig. A. 1(c) does not contain this information, objects have 
been uniformly distributed with respect to the magnitude. 
These facts explain the difference between both figures. 
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where the constants a, b, a, a m in are derived from the 
OGLE 2 data. P(a) represents the probability for having 
the standard deviation a. The constants have been derived 
for each "magnitude bin", for each OGLE fields in SMC 
and LMC. Then average values have been calculated over 
SMC and LMC. 

In our population synthesis code, for a given mag- 
nitude value to, a standard deviation value a m is ran- 
domly determined following the probability law derived 
from OGLE. After that, either the object is rejected (if the 



